Optical Spectroscopy of Tungsten Carbide for Uncertainty Analysis 
in Electron Electric Dipole Moment Search 
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We perform laser induced fluorescence(LIF) spectroscopy on a pulsed supersonic beam of tungsten 
carbide(WC) molecules, which has been proposed as a candidate molecular system for a permanent 
Electric Dipole Moment(EDM) search of the electron in its rovibrational ground state of the X 3 Ai 
state. In particular, [20.6]f2 = 2, v' — 4 <— X 3 Ai,u" = transition at 485nm was used for 
the detection. The hyperfine structure and the fi-doublet of the transition are measured, which 
are essential for estimating the size of the potential systematic uncertainties for electron EDM 
measurement. For further suppression of the systematic uncertainty, an alternative electron EDM 
measurement scheme utilizing the g factor crossing point of the fl-doublet levels is discussed. On 
the other hand, flux and internal temperature of the molecular beam are characterized, which sets 
the limit on the statistical uncertainty of the electron EDM experiment. With the given results, the 
prospect of electron EDM experiment with the X 3 Ai state of WC molecule is discussed. 
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I. INTRODUCTION 

Physics beyond the Standard Model have been devel- 
oped for many years, along with the experimental efforts 
ranging from the Large Hadron Collider to laboratory- 
based tabletop precision measurements [T|. CP viola- 
tion in the Standard Model predicts an electron to have 
a permanent Electric Dipole Moment(EDM) of \d e \ ~ 
1CP 38 e-cm, however, various extensions of the Standard 
Model predict electron EDMs of 10 orders of magnitude 
larger [2j. Therefore, a non-zero measurement of the 
electron EDM becomes a direct observation of physics 
beyond the Standard Model. On the other hand, a null 
measurement of electron EDM with a small enough un- 
certainty would serve as a constrain for theoretical mod- 
els. 

The valence electrons of a suitably chosen diatomic 
molecule have been known to have benefits for the elec- 
tron EDM search Q . The EDM of the valence electron 
gives rise to an energy splitting between opposite spin 
parity states when the molecule is polarized with an lab- 
oratory electric field, £\ a b- This energy splitting is pro- 
portional to the effective electric field experienced by the 
electron, £ c g. There is a big gain of £ e s in heavy po- 
lar molecules compared to atoms Q, which opened up a 
new generation of molecule based electron EDM experi- 
ments In particular, the ytterbium fluoride (YbF) 
system reported a new experimental limit on the electron 
EDM to be \d e \ < 1.05 x 10~ 27 e-cm Q. 

Recently, the A 3 Ai ground state of tungsten car- 
bide (WC) has been proposed as another candidate sys- 
tem for electron EDM measurement Q. This system 
carries two main advantages for measuring the electron 
EDM. The first advantage comes from the large £ c g being 
applied to the valence electrons of the molecule, which is 
due to the heavy tungsten nucleus. The second advan- 
tage is the use of internal comagnetometer H illustrated 
in Fig. [TJ where we can use the closely spaced f2-doublet 
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FIG. 1. Diagram of Zeeman shift cancelation scheme in WC 
electron EDM experiment. The tensor Stark shift /i e ;£iat>, 
Zeeman shift gusB, and the eEDM Stark shift d e £cff are 
shown. 



levels for Zeeman shift cancelation. When g e and gj are 
close enough, we effectively cancel out the Zeeman shift 
terms in Fig. [T]by taking the difference between the en- 
ergy shifts of top and bottom measurements, and only 
leave the Stark shift term, which is proportional to d e . 

On the other hand, there are also uncertainties for two 
of these aspects and therefore careful analysis is required. 
As for the £ e ff, since there is no direct way of measur- 
ing the electric field inside the molecule, the uncertainty 
purely comes from the calculation itself. The calcula- 
tion of £ e ff requires information on electron wavefuction 
at the heavy nucleus of the molecule therefore be- 
ing strongly related to the hyperfine structure. In other 
words, comparing the calculated hyperfine constant with 
the experimental result could be a qualitative test of self 



2 



consistency in theoretical calculation of the £ c g. 

The Zeeman shift cancelation scheme shown in Fig. [1] 
also has room for systematic uncertainty, as the small dif- 
ference in g-factors between the top and bottom doublet 
would result in imperfect cancelation. The magnitude of 
this systematic uncertainty can be written as, 2fj,BBAg, 
where \ib is the Bohr magneton, B is the laboratory mag- 
netic field, and Ag is the difference between g e and gf. 
Petrov [l(| has shown that Ag is closely related to the 
energy splitting between top and bottom levels of the 
f2-doublet. Utilizing this relation, we can calculate Ag 
based on O-doublet constant measurement, and assign 
systematic uncertainty for electron EDM measurement. 

In order to analyze these uncertainties, here we have 
studied the [20.6]ft = 2,v' = 4 <- X 3 A u v" = tran- 
sition of WC by laser induced fluorescence(LIF) spec- 
troscopy on a pulsed supersonic jet of WC molecules. The 
optical properties of WC have been studied both theo- 
retically [ill HI and experimentally [HO]. We had 
~ 100 times higher spectral resolution than the previ- 
ous resonant two-photon ionization(R2PI) spectroscopy 
of the same transition [l3j . The rotational spectrum 
of [20.6]ft = 2,v' = 4 <- X 3 Ai,v v = transition at 
485nm was presented for the low lying R lines(A J = +1). 
The magnetic hyperfine constants and the O-doublet con- 
stants of ground and excited state of the transition are 
reported and analyzed closely related to the systematic 
uncertainties of the electron EDM measurement scheme. 
The flux and internal temperature of the beam are char- 
acterized, which sets the limit on statistical sensitivity of 
the electron EDM experiment. 



II. EXPERIMENTAL METHODS 

We use Smalley type pulse supersonic beam tech- 
nique [18| to generate WC beam with low rotational tem- 
perature. The requirement of low rotational temperature 
is not only to simplify the rotational spectrum, but also 
to enhance the statistics for the electron EDM experi- 
ment, which will be discussed in section V. 

Tungsten atoms are ablated from a rod (American 
Elements, 99.9% purity) by the third harmonic of the 
Nd:YAG pulse laser (Quantel), while the solenoid gas 
valve (Parker, general valve series 999) entrains the atoms 
with 350psi of Argon buffer gas pressure. The WC 
molecules were generated by adding a small fraction of 
methane to the buffer gas, which allows for the chem- 
ical reaction W + CH 4 -> WC + 2H 2 to happen. The 
molecules get cooled down through buffer gas collisions, 
resulting low internal temperatures. The turbo pump 
with 1500L/s of pumping speed maintained the operating 
pressures to be, 5 x 10~ 6 Torr inside the vacuum cham- 
ber. Diagram of the experimental apparatus is shown in 
Fig.H 

A tunable CW diode laser(Toptica DL pro) at the 
wavelength of 484nm-487nm range covered all the tran- 
sitions presented in this paper. The probe laser was 
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FIG. 2. Diagram of tungsten carbide beam apparatus. 

focused with an intensity of ~ 100mW/cm 2 at the in- 
tersection point where the laser beam is crossing the 
molecular beam perpendicularly. The laser induced fluo- 
rescing light was collected by a spherical lens into a wa- 
ter cooled Photo Multiplying Tube(Hammamatsu) con- 
nected to the photon counter. A 485 ± 1.5nm bandpass 
filter was installed in front of the PMT to only let the 
[20.6]ft = 2, v' = 4 X 3 A 1: v" = fluorescence to pass 
through. The photon counts were recorded simultane- 
ously as the wavelength meter(High Finesse WSU series), 
which is a Fizeau interferometer with a frequency refer- 
ence, measures the frequency of the probe laser. The 
uncertainty is defined in two different ways, one is in 
terms of absolute frequency, and the other is the relative 
frequency shift. 

The uncertainty of absolute frequency measurement 
depends on the stability of the frequency reference which 
the wavelength meter synchronizes to. The current set- 
tings use a frequency stabilized HeNe laser(SIOS-02 se- 
ries) with a stability of ~ 10MHz, which was claimed 
by the manufacturer. In other words, there could be 
an inconsistency at the same level in absolute frequency 
measurement. We show on Fig. |3l where the absolute 
frequency of 182 W 12 C, R(l) [20.6]fi = 2,v' = 4 <- 
X 3 Ai,v" — transition was measured over a period of 
several months to show this effect. From the plot, we 
saw an la uncertainty of about 20MHz over 4 month pe- 
riod. This is assigned as a systematic error in absolute 
frequency. 

What we are more interested in is the uncertainty of 
relative frequency shift, as our main concern is to mea- 
sure the difference in energy levels coming from various 
interaction Hamiltonians of WC. If we were to measure 
a relative frequency shift between two LIF lines, the un- 
certainty of each line positions explained above would 
become irrelevant, as it gets canceled out when we take 
the difference in frequency measurements. Therefore we 
are only left with the uncertainty that comes from the in- 
terferometer part of the wavelength meter, where we have 
assigned 1.6MHz of statistical uncertainty and 6 x 10~ 4 
of fractional systematic uncertainty for our relative fre- 
quency shift measurements (see the Appendix of 19]). 
All of these uncertainties are listed on table|H Both statis- 
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TABLE I. Experimental sources of error in frequency mea- 
surement. 



Source of Error 



Estimate 



Statistical Error 

Fractional Systematic Error 
in Relative Frequency 

Systematic Error in 
Absolute Frequency 



1.6MHz 



6 x 10" 



20MHz 



618,111,1 



618,111,060 



618,111,020 



Apr 26 



May 24 



Jun 21 
Date 



Aug 16 



[20.6]Q=2,v'=4, J' 



X 3 A,,v"=0, J" 



« W I2 C 



R(l) 



« W I2 C 



R(2) 



F 

r+i/2 

J'- 1/2 
J"- 1/2 

r+1/2 








FIG. 3. The 182 W 12 C, R(l) line position of [20.6] fi = 2,v' = 
4 <— A 3 Ai,w" = transition has been measured over 4 
months period. 



tical and systematic uncertainties are taken into account 
for la errors of our experimental results in the following 
section. 



III. EXPERIMENTAL RESULTS 
A. Measured WC Transitions 

Tungsten has four major isotopes of 182 W, 183 W, 184 W, 
186 W, while Carbon has only one major isotope of 12 C 
with close to 99% abundance. The 180 W isotope has 
a natural abundance of 0.12%, which was not detected 
with our signal to noise level. As a result, four WC iso- 
topes are observed for each rotational transitions. The 
LIF signal was detected with 1 ~ 10 photon counts per 
second, which we have averaged over 20 second of inte- 
gration time per data point to get better signal to noise. 
The background level is dominated by the randomly scat- 
tered probe laser light going into the PMT. The strength 
of the signal varied on daily basis at 30% level, depending 
on the optimization of the WC beam quality. 

We measured the R branch (A J = +1) of [20.6]fi = 
2,v' — 4 «- X 3 Ai,v" = transition at 485nm, where 



FIG. 4. R(l) and R(2) lines of [20.6]fi = 2,«' = 4 «- 
X 3 Ai,v" = transition was measured for 183 W 12 C. The 
photon count rates are averaged over 20 second of integra- 
tion time, lcr error bars are shown for the photon count rates 
and the red solid line corresponds to the least square Gaussian 
fit. 



the results of 183 W 12 C in Fig. [J] shows the Hyperfine 
structure and the results of 184 W 12 C in Fig. [5] shows the 
fi-doublet structure. We observed the same fi-doublet 
structure in 182 W 12 C and 186 W 12 C isotopes as well, 
which will be shown in the Appendix section. 

Each of the individual line had a Doppler Broadened 
linewidth of ~ 10MHz. Least square Gaussian fit was 
used to extract the center frequency of the line. The fit- 
ting error for the center frequency was only in the order 
of sub MHz, therefore, the uncertainties from the wave- 
length meter discussed in the previous section were the 
dominating sources of error. The list of our measured 
line positions are shown with proper error assignments 
and compared with the previous results of ref. [i~3| on 
Table [Til The center of gravity position is shown for the 
183 W 12 C isotope, as it has a hyperfine structure. For the 
fi-doublet structure observed in R(4) and R(5) transi- 
tions, the center frequency between the double peak is 
listed on the table. Sickafoose's data shows the fitted 
line positions with the residuals from the fit given in the 
parentheses in the units of last significant digit. 



[20.6]Q=2, v'=4, J' 



X 3 A ,v"=0, J" 
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FIG. 5. R(l) ~ R(5) lines of [20.6]f3 = 2, 



X 3 Ai, 



= 4 <- 

transition for 184 W 12 C are shown. The photon 
count rates are averaged over 20 second of integration time. 
la error bars are shown for the photon count rates and the 
red solid line corresponds to the least square Gaussian fit. 



The 183 W 12 C isotope has a non-zero nuclear spin, 
which gives rise to hyperfine structure. Only the mag- 
netic dipole interaction is present as it has a nuclear spin 
of 1/2 [20]. From the selection rule, there are three al- 
lowed electric dipole transitions of AF — 0,±1 for each 
R lines. We have measured 5 out of 6 hyperfme transi- 
tions in 183 W 12 C, R(l) and R(2) lines, which are shown 
on Fig.HJ The 183 W 12 C isotope has a relatively low natu- 
ral abundance of 14% compared to other isotopes, there- 
fore, only up to R(2) transition was detected. Also, the 
183 W 12 C, b transition of R(2) line has a relatively small 
Clebsch-Gordan coefficient which made it undetectable 
with our signal to noise. The intensity ratios of these 
transitions are calculated to be [I a : h : I c = 9 : 1 : 5] for 
R(l) and [I a : I b : I c = 20 : 1 : 14] for R(2). 

Figure [5] shows the five lowest R lines of the 184 W 12 C 
isotope for the given transition. Both the X 3 Ai state 
and the [20.6] SI = 2 have nearly degenerate spin states of 
opposite parity called the S7-doublet, where there are two 
allowed electric dipole transitions among them for each 



TABLE II. Absolute frequencies of measured lines in 
[20.6] fi = 2,v' = 4 «- X 3 Ai,u" = transition are com- 
pared with ref. [ll|. Center of gravity position is shown for 
183 W 12 C isotopes. As for the R(4) and R(5) transitions, we 
start to see an fi-doublet structure, where the center frequen- 
cies between the peaks are listed. The 1 a errors are shown 
inside the parentheses in the order of first parenthesis with 
the statistical uncertainty, and the second parenthesis with 
the systematic uncertainty in absolute frequencies. Only the 
1 a fitting error was provided from ref. [13J 



Measured Line 


This Work 


Ref. [13] 




(MHz) 


(MHz) 



182 W 12 


c, 


R(l) 


618,110,996(1.6) 


(20) 


618,110,980(690) 


182 W 12 


c, 


R(2) 


618,132,466(1.6) 


(20) 


618,132,416(540) 


182 W 12 


c, 


R(3) 


618,150,892(1.6) 


(20) 


618,151,362(-180) 


182 W 12 


c, 


R(4) 


618,166,298(1.6) 


(20) 


618,166,802(-480) 


182 W 12 


c, 


R(5) 


618,178,709(1.6) 


(20) 


618,178,613(-270) 


183 W 12 


c, 


R(l) 


618,103,338(1.6) 


(20) 


618,102,617(-30) 


183 W 12 


c, 


R(2) 


618,124,778(1.6) 


(20) 


618,123,992(-30) 


184 W 12 


c, 


R(l) 


618,093,845(1.6) 


(20) 


618,093,443(540) 


184 W 12 


c, 


R(2) 


618,115,308(1.6) 


(20) 


618,115,208(180) 


184 W 12 


c, 


R(3) 


618,133,727(1.6) 


(20) 


618,133,345(390) 


184 W 12 


c, 


R(4) 


618,149,130(1.6) 


(20) 


618,149,624(-570) 


184 W 12 


c, 


R(5) 


618,161,507(1.6) 


(20) 


618,162,095(-780) 


186 W 12 


c, 


R(l) 


618,077,509(1.6) 


(20) 


618,076,625(660) 


186 W 12 


c, 


R(2) 


618,098,953(1.6) 


(20) 


618,098,509(240) 


186 W 12 


c, 


R(3) 


618,117,361(1.6) 


(20) 


618,117,186(-30) 


186 W 12 


c, 


R(4) 


618,132,760(1.6) 


(20) 


618,132,805(-210) 


186 W 12 


c, 


R(5) 


618,145,138(1.6) 


(20) 


618,145,397(-420) 



R transitions. Following the notations of Brown [21|, 
we label them e/f ■<-> e/f. Due to the J dependence 
of the doublet interaction, HdoiMet °^ J(J + 1), we only 
see the doublet structure when the interaction term be- 
comes larger than our Doppler limited linewidth. As a 
result, we saw the doublet structure only at higher J line 
transitions, which are R(4) and R(5) transitions. 



B. Hyperfme Constant 



As discussed previously, the hyperfine structure of 
183 W 12 C isotope is only cause by magnetic dipole inter- 
action, which can be written as, 



2J(J+1) ' V ; 

where h is the magnetic hyperfine constant. The hyper- 
fine frequency splitting within the electronic state be- 
comes, 
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TABLE III. Relative frequency shifts of measured hyperfine 
splittings of 183 W 12 C in [20.6]fi = 2, v' = 4 <- X 3 A 1 ,v" = 
transition. The 1 a errors are shown inside the parentheses in 
the order of first parenthesis with the statistical uncertainty, 
and the second parenthesis with the systematic uncertainty 
in relative frequencies. The residuals from the fit are shown 
on the third column. 



Measured Line 




Splitting 


Fit Residual 






(MHz) 


(MHz) 


183 W 12 C, R(l), 


|o-6| 


528(1.6)(0.3) 


-2.3 


183 W 12 C, R(l), 


\b-c\ 


882(1.6)(0.5) 


-1.5 


183 W 12 C, R(2), 


\a — c 


849(1.6)(0.5) 


-9.8 



It is straightforward to extract the magnetic hyperfine 
constants of ground and excited state of the transition 
using equation [2j with given hyperfine frequency split- 
tings among the peaks of Fig. 0] First we had to assign 
each measured frequency splittings to the correct hyper- 
fine transition (i.e. whether the splitting corresponds to 
\a — b\, \b — c|, or |c — a\ transition). We have used the 
least square fitting method for all combinations of as- 
signments, and picked the assignment which gave us the 
smallest residuals for the fit. Using the wrong assignment 
made a significant increase in the residuals, which allowed 
us to reject it. The result of three independent frequency 
splittings with the correct assignments are shown on ta- 
ble lllll Based on this assignment, we report the measured 
hyperfine constants to be ft.3Ai,u"=o = — 1171(4)MHz 
and ft.[20.6]f2=2,u'=4 — 1258(6)MHz with la fitting errors 
shown inside the parentheses. Detailed comparison with 
the calculated hyperfine constant will be shown in section 
IV. 



C. fi-doublet Constant 

Both X 3 Ai state and [20.6]fi = 2 state of WC have 
parity eigenstates in each of its J rotational levels, which 
are known as the 57-doublet. These doublet structures 
have been studied by Brown (2lj . where they have pro- 
posed a labeling convention of |e) levels for parity +(— 1) J 
and \f) levels for parity — (— 1) J . The coupling between 
the rotational and electronic motion lifts the degeneracy 
of the parity eigenstates, causingenergy splittings for 
states with 0,^1. Following ref. [21j . we write down the 
Hamiltonian for both states as, 

fl3 Al =±03 Al J(J+l), (3) 

H [20 . 6 ]n=2 = ±5 [20 . 6]n=2 (J - 1) J( J + 1)(J + 2), (4) 



TABLE IV. The ^-doublet constants of 3 Ai and [20.6]fi = 2 
states are shown and compared. Ref. [l7| shows the f2-doublet 
constants of 3 Ai for the 184 W 12 C. The 1 a fitting error is 
shown inside the parentheses for both experimental values. 



fi-doublet Constant 


This Work 
(kHz) 


Ref. [17] 
(kHz) 


03 Al 


418(18) 


400(13) 


O[20.6]n=2 


< 1 





where 03 Al O[2o.6]fi=2 are the O-doublet constants. With 
this, the doublet frequency splitting shown on Fig. [5] can 
be fitted using, 

<5z/ = ±25 3Al J(J+l)±25 [20 . 6] a =2 (J-l)J(J+l)(J + 2). 

(5) 

where the ± signs indicate 4 different possibilities of 
{+, +}, {+, -}, {-,+}, and {-, -} for the fit. We note 
that {—,+}, { — , — } cases would give the same fit result 
as the { + , +}, {+, — } cases, only with the opposite signs 
for the constants. As we only care about the magnitude 
of the doublet constants, this reduces the possibilities to 
2 cases. The least squared fit for two different cases of 
{+, +} and {+, — } both showed small residuals, and also 
revealed that 2o3 Al J(J+ 1) term is much larger than the 
26[20.6]Q=2(y — 1)J(J + + 2) term. In other words, 
the first term was the main contribution to the observed 
frequency splitting, that made the second term to be ir- 
relevant whether it was added or subtracted from the first 
term. 

Fitting the fi-doublet frequency splittings observed in 
R(4) and R(5) lines of 182 W 12 C, 184 W 12 C, 186 W 12 C iso- 
topes gave the ground state f2-doublet constant 03 Al , and 
the upper bound of the excited state fi-doublct constant 
O[20.6]n=2 ; which are shown on table HVl Details of the fit- 
ting procedure, and the LIF measurements of 182 W 12 C, 
186 W 12 C isotopes are shown in the Appendix. The iso- 
tope dependent doublet constants are all within the error 
bars of our extracted constants. The measured fi-doublet 
constant 03 Al is used for the estimation of the difference 
in g-factors between the doublet levels, which will be dis- 
cussed in section IV. 



IV. SYSTEMATIC UNCERTAINTIES OF WC 
ELECTRON EDM EXPERIMENT 

A. Systematic Uncertainty from £ e R 

In the previous section, we have reported on the hy- 
perfine constants of both ground 3 Ai,w" = state and 
excited [20.6]S1 = 2, v' = 4 state. The part that we 
are interested in is the ground state hyperfine constant 
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' l3 A 1 ,t)"=0j as it contains the information of the WC elec- 
tron wavefunction near the tungsten nucleus. As the £ c g 
field comes from the relativistic effect being applied to 
the WC valence electrons near the heavy nucleus @ , the 
calculation of £ e ff in a specific electronic state is closely 
linked to the calculation of hyperfine constant in the same 
state. Therefore, the comparison between the calculated 
WC hyperfine constant /i3 Al v i'—q and the experimental 
results could serve as a qualitative test of self consistency 
in the calculation of £ e g. 

To perform ab initio calculation of the WC 
molecule, a generalized relativistic effective core poten- 
tial (GRECP) [22|, [H| for tungsten atom was generated. 
As a result, 60 core electrons (Is— 4/) of W were ex- 
plicitly excluded from the correlation treatment. Basis 
set for W was constructed using the generalized corre- 
lated scheme [24[ and it consists of 7s-, 9p-, 6d-, 4/- 
and 2g— type contracted gaussians; such a basis is usu- 
ally written as [7,9,6,4,2]. The aug-cc-pVTZ basis set [H| 
reduced to [5,4,2,1] was used for carbon. 

To evaluate £ ff, one needs to compute the following 
parameter of the PT-odd molecular Hamiltonian (dis- 
cussed in Refs. H S& |27| ) : 



W d 



1 



(6) 



where is the wave function of the considered state 3 Ai, 
and Q is the projection of total electronic momentum on 
the molecular axis directed from W to C, 



Hd — 2cL 




ctE 



(7) 



E is the inner molecular electric field, and <r are the Pauli 
matrices. In these designations, £ c g = Wd|f2|. 



The hyperfine constant An for the Ai state is [28 1 



4i 



OLi X Ti 



l*3 Al ) (8) 



where ^ w = 0.11778471/Lt N [29J. 

To calculate matrix elements [6] and we used a two- 
step scheme, where one performs correlation calculation 
for the valence and outer core electrons with the GRECP, 
followed by the nonvariational one-center restoration of 
the wave function at the inner core region of W (see [27j 
for more details). Twenty electron correlation calcula- 
tion was performed using the spin-orbit direct multiref- 
erence configuration interaction (SODCI) approach [30| . 
accounting for the spin-orbit selection procedure [3l[ . As 
a basis set of one-electron functions for the SODCI cal- 
culation, we used the eigenvectors (natural orbitals) of 
some one-electron density matrix calculated at the scalar- 
relativistic coupled-clusters level with sin gle and double 
cluster amplitudes using the CFOUR code [33 | . 

£ e s and Am were calculated at R(W-C)=3.2248 a.u., 
which is close to the equilibrium distance. The calcu- 
lated A\\ value is ~ —1192MHz and the £ e ff value is 



TABLE V. The hyperfine constants of 183 W 12 C in [20.6] = 



2 V 



X d Ai, 



transition is shown and compared. 



The 1 a fitting error is shown inside the parentheses for both 
experimental values. The calculated value is believed to be 
correct within 10% level. 



Hyperfine Constant 


This Work 


Ref. [15] 


Calculated 




(MHz) 


(MHz) 


(MHz) 


h,3 Al 


-1171(4) 


-1363(17) 


-1192(10%) 


ft[20.6]n=2 


1258(6) 







~ — 36GV/cm. The influence of interaction with the low- 
lying electronic state 3 A2 on the hyperfine structure will 
be discuss in the following subsection along with the Q- 
doublet analysis. 

We report the measured hyperfine constants of 
h3A u v"=o and h [20 .6]n=2,v>=4 in Table El For the hyper- 
fine constant of 3 Ai, v" = state, we compare our result 
with the previous results of Wang's [la ], which had a 
slightly larger uncertainty then our case. As seen from 
Table [V] our result of had good agreement with the 
theoretical calculation, however, disagreed with ref. [l5[ 
in 10% level. The disagreement was relatively large con- 
sidering the uncertainties assigned on each experimental 
hyperfine constants. On the other hand, all three re- 
sults of /i3A x constant agrees within 10% level. These 
results give information on the ground state electronic 
wave function at the nucleus of the molecule, which is 
strongly related to the effective electric field calculation 
of the valence electrons in WC. Therefore, we believe our 
effective electric field calculations are self consistent at 
least in 10% level, which is given by the hyperfine con- 
stant comparison. 



B. Systematic Uncertainty from Ag 

The biggest systematic uncertainty of the electron 
EDM measurement with YbF molecules, which holds the 
current experimental limit, came from the imperfect £i a b 
reversal combined with the r.f. phase detuning How- 
ever, this will not be present in our case. Instead of 
reversing the £i a b for the Zeeman shift cancelation, our 
measurement scheme shown on Fig. [T] uses the fi-doublet 
levels of 3 Ai state Q. Our systematic uncertainty would 
come from a different source, which is the small differ- 
ence in g-factors between the top and bottom doublet 
levels. The upper bound of systematic uncertainty in 
our electron EDM measurement scheme can be written 
as 2/j.BBAg. 

As mentioned previously in the introduction section, 
Ag is closely related to the energy splitting between the 
top and bottom levels of the f2-doublet [T(|. We use the 
constant 03& 1 reported in table HVl for the calculation of 
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TABLE VI. Calculated values of HFS as a function of J 
for the 3 A! state of WC. The second(third) and fifth(sixth) 
columns are the results obtained without (with) interaction 
from the 3 A2 state taken into account. 



J 






J 






1 


882.5 


882.5 


16 


72.4 


42.2 


2 


495.5 


492.4 


17 


68.2 


36.2 


3 


347.5 


342.2 


18 


64.6 


30.6 


4 


268.3 


260.8 


19 


61.2 


25.5 


5 


218.7 


209.2 


20 


58.2 


20.6 


6 


184.6 


173.2 


21 


55.5 


16.1 


7 


159.8 


146.5 


22 


53.1 


11.7 


8 


140.9 


125.6 


23 


50.8 


7.6 


9 


126.0 


108.8 


24 


48.7 


3.7 


10 


113.9 


94.9 


25 


46.8 


0.1 


11 


104.0 


83.1 


26 


45.1 


-3.7 


12 


95.6 


72.9 


27 


43.4 


-7.2 


13 


88.5 


63.9 


28 


41.9 


-10.5 


14 


82.4 


55.9 


29 


40.5 


-13.8 


15 


77.1 


48.7 


30 


39.1 


-17.0 



the g-factors of top and bottom levels of the fi-doublet. 
We need to consider the influence of interaction with 
the low- lying electronic state 3 A2 on the hyperfine struc- 
ture and g-factors of the 3 Ai state of WC, where rough 
ten-electron calculations were performed (core states of 
5s 2 5p 6 of tungsten and Is 2 of fluorine were excluded from 
the correlation treatment). The calculated off-diagonal 
electronic matrix elements arc, 



A/2 = fl'<* 3Al |Ji|*3 A3 ) =0.8 cm" 1 , (9) 
T^'E^T^) l* 3 ^> = 2742 MHz,(10) 
G± = {*3 Al \J- + Aa ) =3.1(11) 

The required diagonal electronic matrix elements and 
excitation energies are taken from experiment: B' = 
0.509 cm" 1 , A\i = _im GHz ( this work ). G\\ = 0-022 
[lf|, D — 1.53 a.u. 15], Ez A2 - Ez Al = 1194 cm" 1 0- 

In table IVI1 the hyperfine splittings (HFS) calculated 
from E(F = J - l/2)-E(F = J + 1/2) are given for 
each rotational levels of the 3 Ai state of WC. The hy- 
perfine structure results were obtained by numerical di- 
agonalization of the Hamiltonian, which was described 
in ref.[lC§. Interaction with the 3 A2 state equally influ- 
ences the properties of the |e) and |/) states, whereas 
the interactions with + and 0~ states would lead to the 
different properties between the |e) and |/) states. How- 
ever, the influence of the latter interactions on HFS is 
not considered in the present study due to the complex- 
ity of the WC spectrum. The experimental value of the 
51-doubling allows us to write, 



B ' 2 E 



<*3 A .|^|*„ +) 



E 3 , 



-E 



Er 



(*3 Al \Ji\* n 



0+ 

= ±0.4 MHz. 



E 3/ 



-E_ „_ 



(12) 



Due to a large number of electronic states in the WC 
spectrum, it is hard to reproduce this value in ab ini- 
tio calculations. However, we have estimated from our 
calculations that the matrix elements ($34 1 l^ll^g*) is 
much smaller than ($3^j | J+\ty n Q±) by absolute value. 
Therefore, we can write, 



J2 2^(^3 Al |J^|W n o + )(W3 Al |J^+^|^ n0+ ) _ (13) 



E 



Es Al - E nQ + 
2B'(^ Al \Jl\^ m0 -)(^ Al \Jl + S%\9 m0 -) 



Es Al — E m0 - 



±5 • 10" 5 . 



The ± sign in eqs. (|T2"j) and (fTB")) comes from the fact 
that it has not been experimentally determined whether 
it is |e) or |/) state that belongs to the lower fi-doublct 
level. Taking this into account, the difference of the g- 
factors can be written as \g e — gf \ = 5 • 10 -5 • J(J + 1). 
It can also be shown that interaction with the 3 A2 state 
leads to the same J— dependence of the g-factors for |e) 
and |/) levels: g e(f) = 0.022 + AG±/2(E 3A2 - Es Al ) ■ 
(J + 2)(J-1). _ 

Smaller the difference between g e and gj is, smaller 
the systematics would be coming from spurious magnetic 
fields. This difference depends on the electric field, which 
are shown on Fig. [S] Let us analyze the behavior of g- 
factors in this plot. At first glance, one expects that 
the difference between g e and <?/ could be made zero 
by increasing the electric field, since the external elec- 
tric field mixes |e) and |/) levels. However, this only 
happens for |J = 1,F = 1/2, \M F \ = 1/2) hyperfine lev- 
els of 183 W 12 C, which will be explained in subsection C. 
Without a proper perturbation term from the hyperfine 
interaction, the difference between g-factors would even- 
tually diverge as the electric field increases, which has 
been observed from PbO molecules [33|. This behavior 
was explained by M.G. Kozlov(see acknowledgements in 

Hi). 

Applying the same analysis to the case of WC 
molecules, when J = 2 level is mixed in to J = 1 by 
the external electric field, the g-factor of the lower VI- 
doublet level will have a slight decrease at very low elec- 
tric field, and continuously increase as the field increases 
(here we suggest that electronic g-factor is positive). At 
the absence of external electric field, it was found that 
the initial value for g-factor of the lower fi-doublet level 
is larger than the g-factor of the higher one (34 |. Tak- 
ing these two factors into account, one could understand 
the g-factor curves in Fig. [S][a) for isotopes with zero nu- 
clear spin. Without the hyperfine interaction, the differ- 
ence between g-factors would eventually diverge at higher 
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electric fields, however, was not shown in Fig. EJa) due 
to limited plot range for comparison purposes. 

It can also be shown that whenever signs of the sums in 
Eqs. (fT2j) and ([T3|) are the same, the lower f2-doublet level 
will have a larger g-factor (independently from the parity 
of its state). This is true for WC since ( i i?3A 1 \S+\ i S? ri0 ±) 
is much smaller than ($3Aj | J+I^noi) by absolute value. 
Therefore, we arrive at the same conclusion of g-factors 
diverging at high electric field, for WC isotopes with zero 
nuclear spin. 

On Fig. [5] the calculated g-factors for J=l fi-doublet 
levels of 3 Ai state of WC are given as functions of elec- 
tric field. Here, for simplicity we assume that the value of 
g = 0.022 [l(| belongs to a high fi-doublet level of spin- 
less isotopes. Further experimental investigation would 
be required to confirm this. However, even if it were 
the opposite case, it would only shift all the curves in 
the figure slightly down, leaving their relative positions 
unchanged. 

The maximum EDM induced Stark splitting 2£ c g-d e is 
only reached for fully polarized molecule. For the finite 
electric field, the Stark splitting obtained from numerical 
calculations is plotted in Fig. 

As shown on Fig.^a), we can get the Ag = 0.00002 for 
isotope with zero nuclear spin at £\ a b = lOV/cm, which 
is enough laboratory electric field to fully polarize WC 
molecule. Using this estimation of Ag, we get ~ 100/iHz 
of systematic shift with ~ ^G control of the magnetic 
field. This would limit our EDM sensitivity at 5d e « 
10 -29 e-cm level when £ e g ~ — 36GV/cm. 



T 



T 



low H-doublet zero nuclear spin 
high n-doublet zero nuclear spin 
lowfi-doubletF=1/2 jM f |=l/2 
highn-doubletF=1/2 |M |=1/2 



(a) 



_i i i i i i i_ 



4 6 8 10 12 14 16 18 20 

Lab Electric field (V/cm) 



T 



low n-doublet F=3/2 |M f |=1/2 
highfi-doublet F=3/2 |M F |=l/2 
low n-doublet F=3/2 |M f |=3/2 
highn-doublet F=3/2 |M |=3/2 



8 10 12 
Lab Electric field (V/cm) 



C. Further Suppression of Systematic Uncertainty 
from Ag 

There is an additional interesting feature shown on 
Fig. Oa), where the g-factors cross at £i a \> = 2 V/cm 
for | J = 1,F = 1/2,\M F \ = 1/2) hypcrfinc levels 
of 183 W 12 C isotope due to hyperfine interaction. For 
| J = 1,F = 1/2,\M F \ = 1/2) and |J = 1,F = 
3/2, \M F \ = 1/2) states, the perturbing state is the near- 
est hyperfine energy level, which is simply each other 
(| J = 1,F = 3/2, \M F \ = 1/2) is the perturbing state for 
| J = 1,F = 1/2, \M F \ = 1/2) and vice versa). Note that 
for A {{ <0,\J=1,F= 1/2, \M F \ = 1/2) level is higher 
than | J = 1,F = 3/2, \M F \ = 1/2) level. Therefore, the 
energy denominator for \J = 1,F = 1/2 } \M F \ = 1/2) 
levels in the perturbation theory will have the opposite 
sign compared to the molecule without hyperfine struc- 
ture. Accordingly, if the mixing by the external elec- 
tric field of the | J = 1,F = 3/2, \M F \ = 1/2) level 
is taken into account, the lower Sl-doublet level of the 
| J = 1,F = 1/2, \M F \ = 1/2) will have smaller g-factor 
with increasing electric field. As the initial value for g- 
factor of the lower il-doublet level is larger than the g- 
factor of the higher J7-doublet level, the corresponding 
curves for g e and gf are crossed. On the other hand, for 
| J = 1,F = 3/2, \M F \ = 1/2) level, the perturbing state 



FIG. 6. (a) Calculated g-factor curves for [J = 1,F = 
l/2,\M F \ = 1/2) hyperfine levels of 183 W 12 C, and J = 
1 rotational levels of WC isotope with zero nuclear spin. 
(b)Calculated g-factor curves for |J = 1,F = 3/2, \Mf\ = 
1/2) hyperfine levels of 183 W 12 C, and |J = 1, F = 
3/2, \M F \ = 3/2) hyperfine levels of 183 W 12 C 



is the higher lying \ J = 1,F = 1/2, \M F \ = 1/2) state. 
Therefore the g-factor curves shown in Fig. [6jb) does not 
crossed for the \ J=1,F= 3/2, \M F \ = 1/2) level. The 
perturbing state of | J = 1, F = 3/2, \M F \ = 3/2) level is 
the hyperfine state of the much higher lying J = 2 level, 
which is the \J = 2,F= 3/2, \M F \ = 3/2) state. As a re- 
sult, the g-factor curves of |J = 1,F= 3/2, \M F \ = 3/2) 
level in Fig. [Hb) resembles the g-factor curves for iso- 
topes with zero nuclear spin shown in Fig. |6ja). 

Utilizing the g-factor crossing point of the | J = 1, F = 
1/2, \M F \ = 1/2) hyperfine levels of 183 W 12 C isotope, 
one could further suppress systematic uncertainty coming 
from Ag. There are two main factors to check before we 
could consider using this crossing point for the electron 
EDM measurement. 

The most important factor to check is, how much of 
£ c g field we would get at the the g-factor crossing point 
with the given £i a b- As the crossing point lies at £j a b = 
2 V/cm, which is not enough laboratory electric field to 
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FIG. 7. Calculated EDM induced Stark splitting between 
±M F levels of the \J = 1,F = 1/2); |J = 1,F = 3/2) and 
J = 1 states. All the curves are normalized by the maximum 
EDM induced Stark splitting, 2£ e g ■ d e . The dot-dash line 
shows the EDM induced Stark splitting between ±Mf levels 
of | J = 1, F = 1/2) state at £i ab = 2 V/cm. 



with the EDM measurement with spinless isotope. How- 
ever, this is a small loss compared to the factor of 1000 
gain in systematic uncertainty, therefore would become 
a useful alternative electron EDM measurement scheme 
with enhanced comagnetometer performance. 



V. STATISTICAL SENSITIVITY OF WC 
ELECTRON EDM EXPERIMENT 

Due to the reasons explained in the previous section, 
we believe the WC electron EDM measurement would be 
limited by the statistics rather than the systematic un- 
certainties at this point. In this section, we analyze the 
beam properties in order to show the achievable level of 
statistical sensitivity limit of the electron EDM measure- 
ment, and also discuss about possible ways to improve 
this limit. 



A. Beam Properties Related to Statistics 



fully polarize the WC molecule, we would not get the 
full £ c g ~ — 36GV/cm. However, as shown on Fig. 
£ lab = 2V/cm applied to the |J = 1,F = 1/2, \M F \ = 
1/2) state would still give 85% of EDM induced Stark 
splitting, compared to the fully polarized case. In other 
words, we only loose 15% of the £ c s field, and in return 
we get effectively zero difference in g-factors, which would 
suppress the systematics even further. 

The second factor to consider is the slope of the g- 
factor curves at the crossing point. As we have lim- 
ited fiab stability, the uncertainty in £i a b would get 
transferred to the uncertainty in g-factors at the cross- 
ing point. By simple calculation of, 2x "£i a b field 
stability" x "the slope of the g-factor curve at the crossing 
point" , we get Ag of 0(5) x 10~ 8 , where the number inside 
the parenthesis shows the uncertainty, with ~mV/cm 
control of the £i a b field. Combining the above two fac- 
tors, we expect the systematic uncertainty of electron 
EDM measurement due to Ag to be 1000 times smaller 
than case explained in previous subsection B. 

There is a minor drawback of this EDM measurement 
scheme in terms of statistics. One is due to relatively low 
natural abundance of 14% for the 183 W isotope, which 
is about half of the most abundant 184 W isotope with 
zero nuclear spin. The other drawback comes from a 
relatively small Clebsch-Gordan coefficient for the hy- 
perfine transition. As discussed previously, the intensity 
ratios of the hyperfine transitions in R(l) line follows 
[I a : I b : I c = 9 : 1 : 5], and the \J = l,F = 1/2) 
EDM state can only couple to |J = 2,F = 3/2) excited 
state, which corresponds to the c transition, the second 
strongest. With these two factors, the electron EDM 
measurement scheme using the g-factor crossing points 
of | J = 1,F = 1/2, \M F \ = 1/2) hyperfine levels of 183 W 
isotope would suffer factor of 2 loss in statistics compared 



The spectrum of WC molecule gives multiple infor- 
mation about the beam properties. The height of each 
peaks in Fig. 0] and [5] represent the molecular flux at 
a given electronic, vibrational, and rotational quantum 
state. Conversely, knowing the quantum states and the 
scattering rate of the transition, we can estimate the 
density of molecules within the detection volume. The 
Doppler broadened linewidth can be used for calculation 
of the beam divergence, and the axial velocity can be 
calculated from the time delay between the triggering of 
the gas valve and the PMT gate. 

As the proposed electron EDM measurement scheme 
uses the WC molecules in their rotational and vibrational 
ground state of the A" 3 Ai state [8j, estimation of frac- 
tional rotational ground state population becomes im- 
portant. Intuitively, more molecules would be in their 
rovibrational ground state as the internal temperature of 
the beam cools down by the buffer gas collision. From 
statistical physics, the probability of rotational ground 
state population within the X 3 Ai,v" = state at tem- 
perature T can be calculated from the following equation, 



P(J",T) = 



(2 J" + l)Exp[-^[BJ"(J" +1 



- (2 J + l)Exp[-fe[BJ(J + 1)] 



(14) 

where B is the ground state rotational constant given 
by ref. 13[ . Figure [5] shows the fractional rovibrational 
ground state (v' = 0, J' = 1) population at different 
temperatures ranging from 0.3K to 3000K. From Fig. 0] 
and [5j we observed a general trend of R(l) lines being 
stronger than the R(2) lines. This condition combined 
with equation [14] gives an upper rotational temperature 
limit of T rot < 5.5k. Therefore, we expect minimum of 
35% of WC molecules to be at their rovibrational ground 
state when T rot — 5.5k. For better estimation of the 
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FIG. 8. Fractional ro-vibrational ground state (v = 0, J = 1) 
population at different temperatures ranging from 0.3K to 
3000K. The fraction converges to 1 as the temperature goes 
below IK. 



rotational temperature, accurate ratios among multiple 
rotational lines are required. 

All the parameters mentioned above are listed on ta- 
ble IVIII We separate the directly measured quantities 
from the calculated quantities. In particular, the Frank 
Condon factor has been estimated based on dispersed flu- 
orescence spectroscopy of the same transition [l3| . The 
number of WC molecules within the detection volume has 
been estimated based on the geometric efficiency, quan- 
tum efficiency, scattering rate of the transition, and the 
Frank Condon factor. These parameters will be used for 
calculation of the statistical sensitivity limit of electron 
EDM experiment in the following subsection. 



TABLE VII. The list of measured and calculated parameters 
related to the WC molecular beam. 



Parameter 



Measured 



Calculated 



Beam Divergence 0.00041sr 

Beam Axial Velocity 650m/s 

Distance Between 75cm 
Source and Probe 

Geometric Photon 
Collection Efficiency 

Quantum Efficiency of 
Photo Multiplying Tube 

Franck Condon Factor 

of v' = 4 4— v" =0 excitation 

and v' = 4 — >• «" =0 decay 

WC Molecules 

within detection volume 

T ro t < 5.5K 

Fractional Ro-vibrational 
Ground State 



0.003 



0.1 



0.0001 



10 7 WC/pulse 



> 35% 



TABLE VIII. The list of measured and calculated parameters 
related to electron EDM experiment. 



Parameter 



Measured 



Calculated 



B. Current Level of Statistical Uncertainty and 
Possible Improvements 

The EDM of an electron, d e , inside a fully polar- 
ized WC molecule would produce a frequency shift of 
AfEDM = 2d e £ c ff/h. The statistical uncertainty of the 
EDM shift is defined by the Fourier-limited frequency 
resolution of Az4 ta t = 1/27Tt\/ N x T, where r is the in- 
terrogation time and N is the rate of measurement, and 
T is the duration of measurement. Therefore, we can 
write down the expression for statistical uncertainty in 
the measurement of d e as, 

6d e = k / , (15) 

4ttt£ cS VN x T 

From eq. \TE\ longer interrogation time, larger effective 
electric field, and more number of measurements result 
in smaller statistical uncertainty of d e . Going back to 
Table IVIII one can make the connections between the 
beam properties and the parameters in cq. 1151 Slower 
beam axial velocity would give longer interrogation time 
at fixed distance. Lower beam divergence would increase 



Interrogation Time (r) 1.25ms 

Photon Count Rate (N) 5photons/s 

Effective Electric Field (£ e ft) -36GV/cm 

Projected 8d e , st at 5 x HT 27 ^^ 



the molecular density within the LIF probe laser vol- 
ume, which would result in higher photon scattering rate. 
Colder rotational temperature would also enhance the 
photon scattering rate as more molecules would be in 
their rovibrational ground state. 

The list of measured and calculated parameters in 
eq. [15] are shown on Table IVIIII The current status of 
the experiment gives a projected statistical uncertainty 
of Sd e = 5 x 10 _27 e— cm I y/Day. From our systematic un- 
certainty analysis, we believe this would be the dominant 
term that limits our electron EDM measurement sensi- 
tivity, which is slightly above the current experimental 
limit 0]. 

One possible improvement on this limit could come 
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from extension of the WC beam line to increase the in- 
terrogation time. Due to the EDM sensitive X 3 Ai state 
being in the ground state of WC, our interrogation time 
is only limited by the time of flight of the molecules, not 
by the lifetime of the state. As seen on Table I VIII our 
beam has very low divergence, therefore extending the 
beam line of the molecules could be an option for us in 
order to increase the time of flight. From a conservative 
estimate, we believe factor of 2 ~ 3 of beam line exten- 
sion is plausible without the lost of fluorescence rate. 

Another prospect of improvement for statistics comes 
from the Franck Condon factor. From the dispersed flu- 
orescence spectroscopy of ref. (35|, it was shown that 
v' = 4 ■<— v" — vibrational excitation of [20.6] 57 = 
2 <— X 3 Ai transition happens at ~ 100 times lower rate 
than v' = 1 <— v" = vibrational excitation of the same 
electronic transition. Therefore, probing a new stronger 
transition of [20.6]£7 = 2,v' = 1 <- X 3 A 1 ,v" = at a 
known wavelength of 545nm (l3j |. could give us a factor 
of 100 gain in N. Combining the gains from the vacuum 
system and the laser system, the statistical uncertainty 
of d e is expected to reach 10 _28 e — cm/ y/Day level in the 
near future. 



VI. CONCLUSION 

In this work, we have investigated the electron EDM 
measurement scheme using the X 3 Ai ground state of 
WC molecules, where the focus was on identifying and 
estimating the main systematic and statistical uncer- 
tainties. Well collimated and internally cold beam of 
WC was developed for electron EDM search in X 3 Ai 
of the molecule. We have shown high resolution LIF 
spectroscopy of [20.6]£7 = 2,v' = 4 <- X 3 A 1 ,v" = 
transition to extract molecular constants relevant to the 
electron EDM measurement scheme. 

The ground state hyperfine constant, /i3Ai,«"=0) was 
measured and compared with theoretical calculations to 
show self consistency in effective electric field calculations 
at 10% level. The ground state 51-doublet constant, , 
was measured to give the g-factor curves for top and bot- 
tom doublet levels, and also allowed us to estimate the 
systematic uncertainty coming from imperfect Zeeman 
shift cancelation in our EDM measurement scheme. The 
result leads to systematics of Sd e ps 10~ 29 e — cm level 
with ~ /iG control of the magnetic field. Also an al- 
ternative EDM measurement scheme using the g-factor 
crossing point of hyperfine levels of 183 W 12 C isotope was 
shown, where the systematics could be suppressed even 
further. This crossing point of g-factors only happens 
with molecular isotope with non-zero nuclear spin. 

As for the statistical uncertainty, several beam prop- 
erties such as the flux, velocity, divergence, and rota- 
tional temperature were characterized. These proper- 
ties were linked to the Fourier-limited frequency resolu- 
tion of the EDM induced Stark shift measurement. The 
current statistical sensitivity limit was calculated to be 



TABLE IX. Measured ^-doublet splittings of 182 W 12 C, 
184 W 12 C, 186 W 12 C isotopes in R(4) and R(5) lines of 
[20.6] n = 2,v' = 4 «- X 3 Ai,«" = transition. The 1 a 
errors are shown inside the parentheses in the order of first 
parenthesis with the statistical uncertainty, and the second 
parenthesis with the systematic uncertainty in relative fre- 
quencies. The residuals from the fits are shown on the third 
column for the {+, +} case of eq. [5] and the fourth column 
for the {+, — } case. 



Measured Line 


Splitting Fit Residual 


Fit Residual 






for {+, +} 


for {+, -} 




(MHz) 


(MHz) 


(MHz) 


182 W 12 C, R(4) 


16(1.6)(0.01) 


-1.4 


-1.7 


182 W 12 C, R(5) 


28(1.6)(0.02) 


1.3 


1.5 


184 W 12 C, R(4) 


18(1.6)(0.01) 


0.7 


0.4 


184 W 12 C, R(5) 


27(1.6)(0.02) 


0.3 


0.5 


186 W 12 C, R(4) 


18(1.6)(0.01) 


0.7 


0.4 


186 W 12 C, R(5) 


25(1.6)(0.02) 


-1.6 


-1.4 


Sd e = 5 x 10- 


27 e — cm/ \/ Day. 


Possible improvements 



coming from extension of the beam line and the Frank 
Condon factor were discussed. 

We can summarize the main attractive features of WC 
electron EDM experiment as, (i) verified calculation of 
large £ e ff, (h) use of internal comagnetometer with 57- 
doublet structure, (hi) further suppression of systemat- 
ics due to g-factor crossing point, and (iv) unlimited life- 
time of the EDM sensitive state, X 3 A\. The former two 
features are shared by two other molecule based EDM 
experiments, which are the 3 Ai metastable states of tho- 
rium oxide(ThO) @ and hafnium fluoride ion(HfF + ) [361 ] . 
However, the later two features are unique properties of 
the WC system. Based on our analysis, we believe WC 
system has many advantages for electron EDM measure- 
ment including some unique properties showing a lot of 
promise for future generation experiments. 
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Appendix: O-doublet fitting 

The R branches of [20.6] $7 = 2, v' = 4 <- X 3 A 1 ,v" = 
transition of 182 W 12 C and 186 W 12 C isotopes are shown 
on Fig. [S] and [TO] The 57-doublet splittings observed in 
R(4) and R(5) lines of these two isotopes combined with 
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FIG. 9. R(l) ~ R(5) lines of [20.6]fi = 2,v' 
X 3 A ly v" = transition for 182 W 12 C are shown. 



FIG. 10. 

X 3 A 1 ,v" 



R(l) ~ R(5) lines of [20.6]fi = 2,v' 
= transition for 186 W 12 C are shown. 



= 4 «- 



the results of 184 W 12 C isotope shown in the Fig. [5] are 
listed on the second column of table IIXI 

We couldn't extract any isotope dependence of the 
fi-doublet constant from our data as we did not have 
enough precision on relative frequency shift measurement 
of the doublet splittings. The isotope dependence of 
the fi-doublet constant comes from the difference in re- 
duced mass, which has a fractional difference in ~ 0.1% 
level [13], however, our measured frequency splittings 
have fractional uncertainties in ~ 10% level, which were 
dominated by the fitting errors. This is why we had 
to fit the measured O-doublet splittings from all three 
isotopes using the same fit function of eq. [5] and report 
the isotope independent ground state Sl-doublet constant 
03 / ^ 1 with a relatively large error bar, and only put an 
upper bound of 1kHz for the isotope independent ex- 
cited state fi-doublet constant 6[ 2 o.6]n=2- Nevertheless, 
the extracted 03 Al is consistent with the previous mea- 
surement, which was shown on table HVl 
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